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Pair Production in the Collapse of a Hopf Textnre 
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We consider the collapse of a global “Hopf’ texture and examine the conjecture, disputed in the 
literature, that monopole-antimonopole pairs can be formed in the process. We show that such 
monopole-antimonopole pairs can indeed be nucleated in the course of texture collapse given appro¬ 
priate initial conditions. The subsequent dynamics include the recombination and annihilation of 
the pair in a burst of outgoing scalar radiation. 


I. INTRODUCTION 

An outstanding cosmological problem is explaining 
the homogeneity and isotropy of the universe on large 
scales while allowing for sufficient inhomogeneity and 
anisotropy on small scales to account for the existence 
of galaxies, stars, and, of course, us. Several mechanisms 
for this have been suggested. One of these is the infla¬ 
tionary scenario. An alternative which has also received 
considerable attention is the idea that the production of 
topological defects via phase transitions in the early uni¬ 
verse provided the density fluctuations which eventually 
led to galaxy formation and the like. 

In general, defects result because of a spontaneously 
broken symmetry. If one begins with a set of fields, 
which possess a global symmetry G and interact via a po¬ 
tential V ($) which breaks G down to a subgroup H, the 
vacuum manifold is described by the quotient space G/H. 
Defects will arise in theories for which the vacuum man¬ 
ifold allows nontrivial homotopy groups TTn{G/H) char¬ 
acterized by an integer n. The homotopy groups serve 
to differentiate mappings from the n-dimensional sphere 
S'” into the manifold M. 0. For different values of n, 
one gets different types of defects. Domain walls arise 
for n = 0, global strings for n = 1, monopoles for n = 2, 
and textures for n = 3. 

At very early times, the universe is at a very high tem¬ 
perature so that, in general, the effective potential has 
a unique, symmetric vacuum state which does not allow 
defects. However, as the universe expands, it cools and 
passes through a symmetry breaking phase transition al¬ 
lowing for the possible formation of these defects. The 
vacuum manifold is no longer a single unique state, but 
instead a collection of degenerate vacuum states which 
comprise a nontrivial manifold. The presence of these de¬ 
fects can then serve as gravitational seeds for structure 
formation, and hence the importance of understanding 
their dynamics. 

Here, we will be interested in a particular texture 
model and its dynamics. In accordance with Derrick’s 
theorem, scalar field configurations corresponding to a 
texture are unstable and will tend to shrink or col¬ 
lapse 1^. An interesting conjecture has recently been put 
forth on the collapse of various texture models. Sorn- 
borger et al 1^ suggest that the nature of the various 


collapses can be categorized by whether the particular 
texture model allows more than one type of defect. In 
particular, models for which only is nontrivial (such 
as S'0(4)/S'0(3)) collapse in a similar (and, from their 
results, nearly indistinguishable) manner. On the other 
hand, models for which homotopy groups tts and TT 2 are 
both nontrivial (such as SO{3)/SO{2)) should collapse in 
a qualitatively different manner, with the possible pro¬ 
duction of defects characterized by the other (non-TTa) 
nontrivial homotopy group. These authors provide nu¬ 
merical evidence for their conjecture, having evolved sev¬ 
eral texture models from similar initial conditions and 
examined the collapse process/products. 

The results in would appear to be in agreement 
with interesting experimental evidence in nematic liquid 
crystals. Chuang et al induce phase transitions in a 
nematic liquid crystal described by a SO{3)/SO{2) bro¬ 
ken symmetry and produce textures. These textures, 
they claim, decay via the production of a monopole- 
antimonopole pair. However, later numerical work by 
Rhie and Bennett Q with this model suggests that such 
pair production did not occur. Later still, Luo [| sug¬ 
gests that the reason for the absence of pair production 
in the numerical simulations of Rhie and Bennett is that 
the initial configuration of their simulations does not cor¬ 
respond to that seen in the laboratory experiments of 
Chuang. Instead, he presents different initial data which 
he claims should indeed nucleate monopoles. 

In the following, we analyze the dynamics of collapsing 
Hopf textures, those textures occurring from the break¬ 
ing of SO{3) to SO{2). The results verify aspects of these 
earlier results and clarify the general picture. In the first 
section, we present the model, the evolution equations, 
and discuss the numerical approach. In subsequent sec¬ 
tions, we consider the evolutions of a variety of initial 
data. These include a single monopole and a monopole- 
antimonopole pair which are evolved to demonstrate that 
the code can evolve these objects and to allow the iden¬ 
tification of their presence or absence in later evolutions 
of texture collapse. 
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II. THE EQUATIONS AND NUMERICAL 
APPROACH 

The model we consider is described by the Lagrangian 

= ( 1 ) 

where the potential V (<i>) is 

U($) (2) 

and is a set of three scalar fields. The scalar fields thus 
transform under the global symmetry group S'0(3) with 
the potential breaking this down to SO{2). The resulting 
vacuum manifold is SO{3)/SO{2) = for which the 
second and third homotopy groups are nontrivial. Thus 
this model allows both textures and global monopoles. 
The equations of motion are 

= 47 $“ . (3) 

A trivial rescaling invariance > u allows us to set 
u = 1 , leaving the model with one free, dimensionless 
parameter 7 . 

Following j|] we split the total energy density into three 
parts: the kinetic, gradient and potential energy densities 
defined as 


Pk = \ , 

( 4 ) 

Pa = \ (V$“V4>“), 

( 5 ) 

py =7(4>“$“-1), 

( 6 ) 


where an overdot denotes derivative with respect to time 
and V represents the spatial gradient. The total energy 
density is 


/ cos^ (x) + (2z^/r^ — 1) sin^ (y) \ 

4>(f) = 2 (xzsin^(x)/r 2 +ycos(x)sin(x)/r) 

\ 2 (yz sin^ (y) - x cos (y) sin (y) /r) J 

where r = \J^ y^ ^ is the usual spherical radial 
coordinate and where the boundary conditions on the 
radial function y(r) are y(0) = 0 and y(oo) = tt. We 
discuss the evolution of this initial data as well as the 
evolution of another form of initial data for a texture in 
Sections ^ and 

It is worth noting that for the initial data in Eq. 
the initial total energy density is axisymmetric. For the 
particular case in which y(r) = 2 tan“^ (r), the initial 
energy is actually spherically symmetric. In addition, 
one may expect this axisymmetry to be maintained dur¬ 
ing the evolution. We attempt to take advantage of this 
observation and cast the equations and this initial data 
into a manifestly axisymmetric form. In cylindrical co¬ 
ordinates (p, z,(p), we define new quantities 

= cos p 4)^-I-sin (p (11) 

$“ = sintp — costp <I>^ (12) 


VI. 


which, for the above initial data, can be seen to be man¬ 
ifestly axisymmetric. It is now a simple matter to get 
the equations of motion for Using the equations of 
motion for <I>“ from Eq.(^ we find the equations for 4)^ 
to be 


= 47$+(<i)^$ 
= 47$’($^ 4 ) 


- $+ + $+ 
pz 


+ J_$+ 

p2 ^ p2 


^-v^) (13) 

1 _ _ 2 , 1 _ 

pZ . pZ pZ 

''-u^) (14) 


with 


PT=Pk+Pg + Pv (7) 

This model also has a conserved topological charge. 
The monopole current associated with the charge density 
is 

^abc 4‘^,p 4 ^(^) 

where and Cabc are purely antisymmetric tensors 

of rank 4 and 3 respectively. The conserved monopole 
charge is 

Q = J ko (fx = J Cpx {eijk Cabc <)>y) , (9) 

the integral of the charge density fco ||] • 

Given appropriate initial conditions, these equations 
are now straightforward to integrate in Cartesian coor¬ 
dinates. Some particular initial data that has been sug¬ 
gested for texture collapse is 


+ $+$+ +$ $ . (15) 

Since the evolution will remain axisymmetric, we can dis¬ 
card the terms which have derivatives with respect to p. 

This simplihcation of the problem from 3D to 2D re¬ 
sults in enormous gains from a computational perspective 
and allows us significant improvement in our potential 
resolution of dynamic features in the evolution. However, 
in our code development and tests, we found it very use¬ 
ful to also have a 3D code with which to compare our 2D 
results and results of earlier investigators. 

For those reasons, we have implemented both a 3D 
code and a 2D axisymmetric code using a second order 
Crank-Nicholson finite difference scheme implemented 
with RNPL 0. Out-going radiation conditions based 
on approximate spherical propagation as r —> 00 are im¬ 
posed 

r4>“ = /(t-r), (16) 


2 




for some function /. This condition is very effective in 
limiting reflection from the boundaries. We have verified 
that these codes are fully convergent and that they are 
stable for many crossing times. 


III. SINGLE MONOPOLE 


As stated in the introduction, we want to evolve a 
collapsing Hopf texture and examine the possible nucle- 
ation of monopole-antimonopole pairs. However, in order 
to understand the behavior of any nucleated monopoles 
modeled by these codes, we first examine the evolution of 
initial data explicitly containing a monopole. Using the 
standard hedgehog ansatz, we can represent a monopole 
of unit charge at the origin by 


where 

/(O) = 0 /(oo) = 1. (18) 

To find the static monopole /(r), we substitute the 
ansatz ( 0 ) into Eq.(^) and set all time derivatives to 
zero 

/" + ^-^' = 47/(f-l). (19) 

r \ / 


We then solve to find f{r). However, instead of using 
the exactly static solution, /(r), which is exact only in 
the continuum limit, we note that the static monopole is 
rather well approximated by 

/(r) « tanh (^r). (20) 


Using this as our initial data shows both the energy and 
monopole charge to be well conserved. 

Because the monopole is a topological object, we ex¬ 
pect that any initial data which has the same boundary 
conditions will also have unit topological charge. Hence, 
we can choose the function f{r) different than /(r) but 
still obeying Eq. (HD to set our initial data. This data 
will not be static but will still have a conserved monopole 
charge. As an example, we examine the evolution of ini¬ 
tial data of the form 


/(O 


1 

2 


1 -I- tanh 


r — R 

~T~ 


( 21 ) 


with R and (3 arbitrary constants. 

These simulations have been done with time symmetric 
initial data. However, we have also implemented a boost 
of the monopole by giving it an initial velocity in some 
direction along an axis. The monopole is seen traversing 
the numerical domain. That we can evolve a monopole 
for many crossing times and also boost it while retain¬ 
ing its particle nature, provides further evidence that the 
global monopole is stable to radial perturbations as dis¬ 
cussed in |8- W . 


IV. MONOPOLE-ANTIMONOPOLE PAIR 


Next, we would like to consider evolving a monopole- 
antimonopole pair with our codes. An anti-monopole is 
just a reflection of the monopole ansatz (^. Thus, it 
takes the form 

= , 22 ) 

Strictly speaking, there is not a superposition principle 
for monopoles and anti-monopoles. However, it is rea¬ 
sonable to assume that if the separation between the 
monopole and anti-monopole making up a pair is con¬ 
siderably larger than the core width of the respective 
particles, a pair configuration can be well approximated 
by matching a monopole solution to an anti-monopole 
solution. To achieve this, we translate the monopole up 
and the antimonopole down some distance L to get the 
solution for a monopole-antimonopole pair located on the 
2 axis 

y ) , (23) 

\z\-lJ 

where, in analogy with the previous section, we consider 
two forms for / 


$ (x) = 


Jx'^ +y'^ + (l^l - Lf 


f = tanh 



and 


(24) 



1 -I- tanh 


^Jx^+y^ + {\z\-Lf-R 
/? 


(25) 


At z = 0, the fields though not differentiable, are 
continuous. This non-smoothness gets smoothed quickly 
by the numerical evolution and does not appear to affect 
the evolution of the pair significantly. 

The total charge of this system during the evolution 
should be identically zero. Indeed, we confirm this nu¬ 
merically as the total charge remains essentially zero 
throughout the evolution. Though a necessary check for 
our code, this fact makes it a bit more difficult to follow 
the evolving charges individually. One quantity which we 
do track is the topological charge density. This, of course, 
should be peaked at the locations of the monopoles. How¬ 
ever, we also found it useful to define and calculate a 
quantity which we call the “half-space charge.” It is sim¬ 
ply the topological charge density integrated over positive 
z only 


Qi/2 — 


dx dy dz. = 


X J y J z—0 


p J z—0 


2'Kp dp dz. 
(26) 
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We emphasize that this quantity has no topological 
meaning and is not necessarily conserved. However, as 
with our earlier argument for superposing the monopole- 
antimonopole pair, if we imagine the constituents of the 
pair to be well separated and individually well localized, 
the half-space charge should approximate the charge of 
the monopole (or anti-monopole). 

We then evolve the pair conhguration for L = 15 and 
show the resulting half-space charge in Fig. Note that 
for early times, the average value of this quantity for 
large enough 7 is roughly unity as we would expect for an 
isolated monopole. For smaller 7 , the half-space charge 
becomes zero much more quickly. This behavior would 
appear consistent with the fact that decreasing 7 has the 
effect of increasing the size of the respective monopoles. 
The increased size of the core makes it more difficult for 
us to associate the half-space charge with the true topo¬ 
logical charge and results in more rapid pair annihilation. 


cal invariant. For it to be so, we would necessarily have to 
extend the integral out to infinity which we clearly can¬ 
not do while distinguishing the monopole from the pair. 
As further support that these oscillations are represen¬ 
tative of the dynamics, oscillations of the total potential 
energy are observed to be correlated to the oscillations 
in the half-space charge. 

In Fig. 1^, the results of varying the parameter L in 
the pair configuration are shown. The graph of the half¬ 
space charge shows that the lifetime of the pair is clearly 
proportional to the initial separation. The bottom graph 
displays the position on the z axis of the maximum charge 
density as a function of time. Because of our somewhat 
ad hoc matching procedure, the monopoles are not influ¬ 
enced by each other until a sufficient amount of time has 
passed for a signal to have propagated between them. At 
that point, the monopoles experience an attractive force, 
approach each other and annihilate. 
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FIG. 1. The half-space charge for the evolution of a 
monopole-antimonopole pair. The monopoles are described 
by Eq.(p5|) with initial separation 2L = 30, (3 = 3.5, 
R = 10 on a 257x513 grid with outer boundaries at 
Zmax = —Zmin = Pmax = 40. For Small 7 , the half-space 
charge quickly goes to zero, while for larger 7 a well defined 
annihilation occurs releasing a large amount of bosonic ra¬ 
diation. The charge for the other half-space is simply the 
reflection of the graph about zero (the dashed line). 

We have investigated the oscillations appearing in the 
half-space charge seen in Fig. |^. The oscillations con¬ 
verge with increasing resolution. Further, by moving the 
boundary farther from the monopole pair while keeping 
the same effective resolution, the oscillations remain the 
same. Both these results indicate that the oscillations are 
a real phenomenon and not a numerical artifact. Hence, 
the oscillations appear to be part of the dynamics. That 
this quantity oscillates is not inconsistent with charge 
conservation as this half-space charge is not a topologi- 
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FIG. 2. Evolutions of various separations L for a 
monopole-antimonopole pair. The top graph shows the 
half-space charge. The lifetime of the pair is expectedly de¬ 
pendent on the separation. The bottom graph shows the po¬ 
sition on the z-axis of the maximum charge density up to 
the time of annihilation. Taking this value to represent the 
location of the monopole, the graph shows the monopole ac¬ 
celerated towards annihilation after some time dependent on 
the initial separation. These were run on a 65x129 grid with 
7 = 0.1. Ghanges to R and /3 change these graphs very little 
(here, R = 1.5 and (3 = 1.0). 


V. “SPHERICALLY SYMMETRIC” TEXTURE 


With a good idea of what the signature of a monopole- 
antimonopole pair might be, we now evolve the initial 
data given in Eq. with 


X(f) = 


1 -I- tanh 


r — R 

~T~ 


(27) 
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Though the precise shape of this function does not appear 
to significantly affect the results, the function x(r) used 
here differs from that of [|j because the function shown 
there is not smooth. 

As we stated in Section the total energy density is 
spherically symmetric only for x(r) = 2 tan“^(r). We are 
therefore being a bit careless to refer to this as spherically 
symmetric initial data. In addition, as shown in Q , sym¬ 
metric configurations in the model we are considering are 
in fact unstable to nonspherical collapse. Nonetheless, 
we will continue to refer to the initial data described by 
Eq.(|T^ and Eq.(p7|) as “spherically symmetric” mainly 
to distinguish it from another set of initial data which we 
investigate that has toroidal symmetry. 

On evolving this data, we find that as the texture col¬ 
lapses, the energies evolve as described in [Q. The kinetic 
and gradient energies evolve towards equipartition while 
the potential eneray remains a relatively small fraction 
of the total. Fig. g shows the energy versus time for the 
collapse of a typical “spherically symmetric” texture. 

However, our results in both the 2D axisymmetric and 
the 3D codes do not bear out the claim that a monopole- 
antimonopole pair is formed [^. As one indication, the 
half-space charge of the “spherically symmetric” texture 
as shown in Fig. ^ does not resemble that shown for the 
pair in Fig. |^. Fig. ^ suggests that there is a separation 
in the monopole charge density, but the duration of the 
resulting peaks in the charge density is independent of 7 
and the constant R. Further, the half-space charge does 
not asymptote to ±1 as one might expect were a pair to 
form and move apart. This provides additional evidence 
that no monopoles are forming. 

The half-space charge is also plotted for variations 
of R and [3 in Fig. ||. Again, the figure suggests that 
the “spherically symmetric” texture does not nucleate 
monopole pairs. 

One measure, given in and Q], for determining if a 
pair is present is that the potential energy density should 
be sharply peaked at the location of a monopole or an¬ 
timonopole. We do see some peaking in the potential 
(and total) energy density for a brief period on the axis 
of symmetry, but the peaks do not persist. Indeed, these 
peaks quickly diminish in size and propagate off the grid 
with the remainder o f the outgoing radiation (see further 
discussion in Section VII ). 

Thus we confirm earlier results of Rhie and Bennett Q 
that pair production does not occur with this initial data 
in contradiction to that claimed in IHl. 



Time 


FIG. 3. The energy components for the “spherically sym¬ 
metric” texture as a function of time. Here, 7 = 0.1, R = 10, 
and /? = 3.5 with outer boundary Zmax = |^min| = Pmax = 40. 
This run was on 513x1025 grid. The inset shows the total 
energy as a function of time. The total energy is conserved to 
less than one percent until gradient energy reaches the outer 
boundary. 



FIG. 4. The half-space charge as a function of time. Initial 
data is the spherically symmetric texture with A = 10 and 
/3 = 3.0. Run on a 129x129x129 grid. 
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FIG. 5. Half-space charge for the “spherically symmetric” 
texture for different values of R and /3. The top graph shows 
the results from variations of j3 with R = 13. The bottom 
shows the results of variations of R with (5 = 3.5. Here, 
7 = 0.1 on a 513x1025 grid. 


<i>(a:) = 


sin I 27ra 




cos ( 27ra ^^ 


r-^+a^ 

—a^'^x-\-2azy 


) y'((r2-o2)2+4a2z2)(2,2+j^2 

Sin ( 2 ™ - -{r--a-)y+2a.. 


+“ ) y'((r2-o2)2+4a2z2)(2,2+j^2) j 

(33) 


The evolutions of the toroidally symmetric tex¬ 
ture generically display the creation of a monopole- 
antimonopole pair for all the values of a we attempted. 
Fig. II shows the half-space charge for a typical run. 
These results are strikingly similar to those shown for 
the explicit pair evolved in Fig. |^. The charge begins at 
zero, quickly reaches a value approximately ± 1 , and re¬ 
mains there until collapse. As 7 is decreased, the behav¬ 
ior fundamentally changes in the same manner as that for 
the explicit pair. Hence, this charge separation appears 
intimately connected with the potential as would be ex¬ 
pected for pair production. This connection is in contrast 
with the nearly 7 -independent results of the “spherically 
symmetric” texture in Fig. 


VI. TOROIDALLY SYMMETRIC TEXTURE 

That the “spherically symmetric” texture does not col¬ 
lapse via the production of global monopoles does not ex¬ 
clude the possibility of pair production, merely our choice 
of initial data. Indeed, experiments with liquid crystals 
provide evidence that monopole pairs are in fact formed 
in Hopf texture collapse. For this reason, different ini¬ 
tial data is suggested in [|| as possibly leading to pair 
creation. This initial data is toroidally symmetric. It is 
given by 


$i=cos0(r?) (28) 

$2 + ^ sm9{rj) (29) 

where, in terms of Cartesian coordinates {x,y,z), the 
toroidal coordinates {ri,l3,(p) are defined as 


tanh 77 = 

2a\Jx^ + Ip- 

(30) 


-b -b z + 



2az 


tan/3 = 


(31) 


a;"' -b 2 / + z — 


tan (75 = 

y 

(32) 


The parameter a is the value of the so-called degenerate 
tori in this coordinate system. The function 9{ri) has 
boundary conditions 0 ( 0 ) = 0 (on the axis of symmetry 
and r 00 ) and 0(oo) = tt (the degenerate torus). Thus 
a simple choice for 0 is 9{r]) = Trtanh)?^). Hence, we can 
write the initial data as 
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FIG. 6. The half-space charge as a function 

of time. The initial data is the toroidally symmet¬ 
ric texture with a = 10.01 and with outer boundary 

^-max/min ^max/min ymax/min i40. RuU OU a 

129x129x129 grid. 

Even more convincing evidence for the production of 
a monopole pair resides in actual movies made from the 
evolutions. Looking at either the potential energy den¬ 
sity or total energy density, two peaks representing the 
monopoles occur along the axis of symmetry. The peaks 
leave the center of the grid, travel outward on the axis, 
and eventually switch direction and annihilate. 
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VII. DISCUSSION 

The study of topological defects is important for stud¬ 
ies of structure formation in cosmology, studies of su¬ 
perconductivity, liquid crystals, and superfluidity in con¬ 
densed matter, and studies of regularity in topology. 



0 10 20 30 40 0 10 20 30 40 


time 

FIG. 7. Comparison of the axisymmetric and three dimen¬ 
sional codes using the toroidally symmetric texture as ini¬ 
tial data. The 2D runs (open circles) were on a 513x1025 
grid while the 3D grid (crosses) was 129x129x129. The graph 
shows agreement between the two codes. 

Here, we construct a code in 3D and in axisymmetry 
to model defect dynamics. Tests of the both codes show 
that they converge with increasing resolution and con¬ 
serve energy. Further, the 3D code duplicates the results 
generated with the axisymmetric code as shown in Fig. |^. 

Our code clarihes some of the conditions under which 
a Hopf texture nucleates a monopole-antimonopole pair. 
Specihcally, our results confirm those of Rhie and Ben¬ 
nett 0 and Luo Q that the “spherically symmetric” tex¬ 
ture does not produce monopole pairs. This hnding is 
in contradiction to the conclusion of However, our 
results are in agreement with their conjecture that mod¬ 
els of texture collapse which allow more than one non¬ 
trivial homotopy group will evolve differently from those 
with just non-trivial tt^. That toroidally symmetric ini¬ 
tial data in this model leads to pair production agrees 
with their observation that these models can nucleate 
defects characterized by the non-TTa homotopy group. 

We speculate that the reason monopole production was 
claimed in |^] is that the collapse of a “spherically sym¬ 
metric” texture does produce massive radiation in shells 
which are peaked along the axis. By examining an isosur¬ 
face of the potential energy, apparent balls are seen ra¬ 
diating along the axis of symmetry. However, our evolu¬ 
tions suggest that observing isosurfaces can be deceiving. 
Indeed, by examining the charge density, no localization 


is observed. 

Consider, for example, Figs. and ^ These show the 
charge densities for four snapshots in time for both the 
toroidally symmetric and “spherically symmetric” tex¬ 
tures evolved using the 2D axisymmetric code. Examin¬ 
ing the toroidally symmetric texture hrst, pair nucleation 
is quite dramatic between the first two frames. From 
an initially vanishing charge density, two localized, op¬ 
positely charged regions move outward from the origin. 
Between the last two frames, the monopoles change di¬ 
rection and move inward towards annihilation. 

Contrast these dynamics with that shown in Fig. || for 
the “spherically symmetric” texture. Peaks in the charge 
density appear on the axis near the origin, but at no 
time are localized regions of non-vanishing charge density 
seen moving outward. Instead, the dynamics show only a 
region surrounding the origin in which the charge density 
oscillates. This is consistent with the form of the half¬ 
space charge in Figs. ^ and which shows a steep rise 
and fall in this quantity. 

In addition, we show the potential energy density as¬ 
sociated with these evolutions in Figs. and |^. If 
monopoles exist in the space, then regions of non-zero 
charge density should be accompanied by trapped po¬ 
tential energy. For the toroidally symmetric texture, 
Fig. |l^ demonstrates that the localized regions of non¬ 
zero charge density in Fig. also contain localized po¬ 
tential energy. Note that these localized peaks in the 
potential energy maintain their distinct character until 
they recombine and annihilate. Hence, we conclude that 
this provides strong evidence that these localized regions 
represent monopoles. 

However, for the “spherically symmetric” texture 
shown in Fig. |^, the potential energy remains concen¬ 
trated around the origin while shells of radiation move 
outward. Again, there is peaking in the potential energy 
density, however these concentrations do not maintain 
any distinct form. In particular, at times there are two 
peaks and at others there is only one in the strong field 
region near the origin. Outgoing radiation can be ob¬ 
served (this is seen more clearly in the evolutions) which 
is peaked along the axis, but these are not accompanied 
by a similar peaking in the charge density. Taken to¬ 
gether, we therefore conclude from this evidence that no 
monopoles are nucleated with this initial data. 

As we mentioned earlier, supplementary to these 
plots, we have produced movies of some of these 

evolutions. _ These evolutions can be viewed at 

http://godel.ph.utexas.edu/'ehirsch/hopf.html. 

Both the axisymmetric and 3D codes are written in 
RNPL, an advanced numerical language developed at the 
Center [Q. Development of this language is ongoing, and 
we expect soon to have the ability to take advantage of 
automatic parallelization so that higher resolutions can 
be run. We also expect to adapt these codes to model 
other broken symmetries. 
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FIG. 8. Evolution of the monopole charge density for a 
toroidally symmetric texture (a = 10.01). A monopole pair 
is nucleated by t = 7.5. By t = 27.5, the pair is accelerated 
towards each other about to annihilate. 



FIG. 9. Evolution of the monopole charge density for a 
spherically symmetric texture {R = 5 and /3 = 3.5). The 
charge density near the origin separates such that there are 
positive and negative regions. However, at no time do these 
regions move apart, and instead they oscillate near the origin. 
In addition, no correspondence between regions of non-zero 
charge density and trapped potential energy (as shown in 
Fig. is observed. 
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FIG. 10. Evolution of the potential energy density for the 
toroidally symmetric texture shown in Fig. g. At times t = 7.5 
and t = 17.5, localized regions of trapped potential region 
correspondii^ to the regions of localized charge density as 
shown Fig. H are observed. These regions appear to repre¬ 
sent monopoles. At t = 27.5, the monopoles are about to 
annihilate. 






FIG. 11. Evolution of the potential energy density for the 
spherically symmetric texture shown in Fig. Potential en¬ 
ergy is seen near the origin. Massive radiation is seen ema¬ 
nating from the dynamics at the origin, however, no trapped, 
localized regions which might correspond to monopoles are 
observed. 
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